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Interference effects in quantum dots between different transport channels can lead to a strong 
suppression of conductance, which cuts like a canyon through the common conductance plot [Phys. 

Rev. Lett. 104 , 186804 (2010)]. In the present work, we consider the thermoelectric transport prop¬ 
erties of the canyon of conductance suppression using the second-order von Neumann approach. 

We observe a characteristic signal for the zeros of the thermopower. This demonstrates that ther¬ 
moelectric measurements are an interesting complimentary tool to study complex phenomena for 
transport through confined systems. 

PACS numbers: 72.20.Pa, 73.63.Kv, 73.23.HK 


I. INTRODUCTION 


Thermoelectric effects in nanoscale structures have 
been frequently studied with the aim to improve the ef¬ 
ficiency of devicesIn addition, it has been recently 
demonstrated that thermopower measurements can serve 
as an interesting tool to characterize complex scenar¬ 
ios due to coherences and interactions in nanoscale sys- 
tems.l^Hni Here the intrinsic advantage of thermoelectric 
measurements is that they probe asymmetries around the 
Fermi level and therefore can easily provide information 
about excited states 

Close to a degeneracy of energy levels in a quantum 
dot, interference and correlation effects can play an im¬ 
portant role for transport and lead to pronounced quan¬ 
tum mechanical phenomena. For a spin-degenerate quan¬ 
tum dot (QD) the conductance experiences an enh ance- 
ment for low temperatures due to the Kondo effect 
The thermopower of the single and multiple quantum 
dots in the prese nce o f the Kondo effect was examined 
both theoreticall}®^^^ and experiment ally.On the 
other hand, in the case of two degenerate levels with equal 
spin, at degeneracy, conductance can be suppressed due 
to electron correlation and interference effects.l^^Jf^ Such 
a spin-polarized two-level model was widely used to in¬ 
terpret the phenomena of phase lapses of transmi ssion in 
Aharonov-Bohm interferometer containing a 
Additionally, due to their small size the quan tum dots 
also can exhibit Coulomb blockade effectThe ther¬ 
mopower of the usual Coulomb blockade sequential tun¬ 
neling peaks and the cotunneling signal were addressed 
in Refs. |44H52j 

In Ref. [53] a system of two spin-polarized degenerate 
levels was realized in an InSb nanowire QD, where dif¬ 
ferent ^-factors allow to control level crossings for the 
same spin by a magnetic field. The experiment showed, 
in good agreement with supporting calculations, that the 
suppression cuts as a canyon through the standard con¬ 
ductance plot for different parities of the level couplings. 
A more detailed analysis of the conductance spectrum 
can be found in Ref. (54] A related two-level system was 
optimized for achieving high thermoelectric performance 
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Fig. 1. (Color online) Schematic of the system: A quantum 
dot with single particle levels Ei and Ei> coupled to the leads 
via tunneling barriers. The leads have a temperature differ¬ 
ence AT which can give rise to a current flowing through the 
dot. 


in Ref. [55l In this work we further elaborate on ther¬ 
moelectric properties and focus on the fingerprint of con¬ 
ductance suppression in the thermopower signal. 


The paper is organized as follows. In Section [H] the 
model for the spin-polarized two-level quantum dot is in¬ 
troduced. Results for conductance and thermopower are 
presented in Section [TTT| Here we focus on the zeros of the 
thermopower, which are relatively easy to extract exper¬ 
imentally. We start in Section [Til A| with a discussion of 
the non-interacting case, where results are obtained using 
transmission formalism. Here we show that up to five ze¬ 
ros in the thermopower can be found if the gate voltage is 
varied. Furthermore, we establish the role of temperature 
and level broadening for the existence of multiple zeros. 
The impact of finite QD charging energy is addressed 
in Section HIB[ where calculations are performed by the 
second order von Neumann (2vN) approach. This reveals 
the full scenario for the canyon of conductance suppres¬ 


sion. Concluding remarks are given in Section IV 
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II. MODEL 

The two-level quantum dot system under investigation 
is depicted in Fig. Its energy levels Ei, Ey are as¬ 
sumed to be individually tunable. In practice this is 
achieved by a gate voltage Vg = {Ei + Ey — iiy — iiy)/2e 
(where e < 0 is the electron charge, fiy/R is the chemi¬ 
cal potential in the left and right lead, respectively); and 
by controlling the detuning /S.E = Ei — Ey , which, for 
example, can be physically realised by a magnetic field if 
the two levels have different ^f-factors.l^ 

The total Hamiltonian for the system consisting of the 
QD coupled to two leads via tunneling barriers can be 
written as 

H = ( 1 ) 


linear relation between AT and the bias V for which this 
current vanishes. This defines the Seebeck coefficient (or 
thermopower) S via 


This coefficient can be numerically obtained for a given 
transport model. For the units of the thermopower we 
use So = ky/lel ~ 86.1 /iV/K and for the units of con¬ 
ductance we use Go = ~ 38.7 /iS. Lastly, in linear 

response to applied bias V or temperature difference AT 
the Hamiltonian Eq. 0 can be mapped to generalized 
Anderson model, which was discussed in Ref. |34l This 
mapping for model parameters considered in this paper 
is discussed in Appendix [B] 


The two-level spinless QD Hamiltonian Hy,, in a single 

,. 1 1 . . A. Non-interacting; case T = 0 

particle basis, is ^ 


Hy — E-^d^d-^ T T U djd^^dj/d^/, (2) 

where U is the charging energy due to Coulomb repulsion 
between electrons in the dot. In our calculations and 
corresponding plots we shift the zero value of the gate 
voltage to the particle-hole symmetric point and use 

Vg = \e\vg - U/2. (3) 


For the non-interacting case, when there is no charg¬ 
ing energy 1/ = 0, the thermopower S and conductance 
G can be calculated exactly from the electronic transmis- 
Sion T{E) aPHsHl 


G = e^Lo, 


S 


1 y 

\e\TLo’ 


(8a) 

(8b) 


The lead Hamiltonian Hyy is 


with 


^LR — 

k,£ 


(4) 


- i 

k J —oo 


dE{E - lif 


(_M)r,E),,., 


where k denotes the lead state and = L, R denotes the 
lead. We assume, that the leads are thermalized resulting 
in an occupation function 

^ e{Ek-tie)/kBTe ^ 

with different temperatures Tl/r, for the left and right 
lead, respectively. The tunneling between the dot and 
the lead is governed by the Hamiltonian Hy that reads 


where the limit of vanishing temperature difference 
AT ^ 0 and vanishing bias V ^ 0 between the leads 
was taken. The transmission is calculated using the Car- 
oli formulsP^ 

T{E) = Tv[yG^{E)TRG^{E)], (10) 

where G^^^{E) is the retarded/advanced Green’s func- 
tion of the quantum dot electrons and F^ is the coupling 
strength given by 


+ h-c-, (6) 

k,£ 






^£,ij — 


( 11 ) 


where, for simplicity, the tunneling amplitudes tu are as¬ 
sumed to be the same for all values of k and the coupling 
strengths are defined as T£i{E) = 27r|t^ip 6{Ek—E) « 
27TUF\t£i\‘^. Here we performed the k-sum using the flat 
density of states approximation, i. e., ^ uy dT, 

with i£f denoting the density of states at the Fermi level 
and 2D denoting the bandwidth of the leads. We also 
assume 2D to be the largest energy scale in the problem. 

In a nonequilibrium situation a current / through the 
dot can be generated both by an applied bias V = 
(/iL — /^r)/^ and a temperature difference AT = Tl — Tr 
between the leads. Close to equilibrium this establishes a 


The retarded/advanced Green’s function can be obtained 
from 


G=(Go-i-E)-\ 


G 


-1 

0 


(z-Ei 0 \ 

0 z-Ev)' 


(12a) 

(12b) 


by replacing z = E + yq iox the retarded function G^ 
and z = E — iq for the advanced function G^, where q is 
positive infinitesimal. Here E is the self-energy given by 

E = El + Er, Ef ^ = Tiy, (12c) 
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in the case of infinite bandwidth D +oo. If the fol¬ 
lowing choice of tunneling amplitudes (corresponding to 
the canyon of conductance suppression) is made 


tni = t, tRi = t, tLi' = —at, tRi/ 
we get the coupling strength matrices 


Ft =r 


-a a 


rR = r 


1 “ha 
-Fa a^ 


at, (13) 


(14) 


with r = 27rz/i?|tp, and the transmission function takes a 
simple and intuitive form as a sum of two Breit-Wigner 


resonances* 


T{E) = 




1 


E — E\ tr E — E\f -|- ia^r 


(15) 


Here the minus sign between both terms relates to the 
different parities of the tunnel couplings for both levels 
as defined in Eq. (13). For small E both terms cancel at 
Ey = o?Ei resulting in zero transmission. This provides 
a line of conductance suppression in the {Ei^Ey) plane. 

In order to see whether the canyon of conductance sup¬ 
pression gives qualitatively different results for the ther¬ 
mopower, we will make a comparison to the case when 
all the tunneling amplitudes have the same sign: 

tni = t, tRi = t, tLi' = +at, tjii' = at. (16) 
This choice leads to such coupling strength matrices 


Fl/r = F 


1 “Fa 
“Fa a^ 


(17) 


The transmission function acquires a more complicated 
structure than a sum of two resonances 


T{E) = F^ 


a^{E - El) + {E - El.) 


{E - El ^ ir){E - Ey -F la^T) + a^P 


(18) 

For zero energy F = 0 the above transmission is zero at 
Ey = -a^Ei. 


B. Second order von Neumann (2vN) approach 

For the interacting case , w e use the second order von 
Neumann (2vN) approach,!^ which contains both coher¬ 
ent effects and cotunneling. However, the method fails 
around and below the Kondo temperature where strong 
correlations between the dot and the lead electrons ap¬ 
pear. For our considered spinless system the Kondo tem¬ 
perature can be defined through the mapping to gener- 
alized Anderson model as discussed by Kashcheyevs et 
q^/| 34 | 67 | Numerically, the thermopower is found in an it¬ 
erative procedure by varying the voltage until the current 
vanishes. The governing equations of the method and the 
solution procedure of those equations is presented in Ap¬ 
pendix In the numerical calculations we use kyT ^ F 


and a finite temperature difference AT = Tl — Tr (about 
an order of magnitude lower) in order to achieve conver¬ 
gent and reliable results. Furthermore, we define the av¬ 
erage temperature T = (Tl + Tr)/2 as a reference point. 
All the calculations in Figs. |4|[7| are done using the 2vN 
approach. For the non-interacting case t/ = 0 the 2vN 
appro ach re produces the exact transmission formalism 
resultWe explicitly checked that for present cal¬ 
culations the values ky/ST = V = O.IF give the linear 
conductance G and the linear thermopower S on the scale 
of Figs. izl within the 2vN approach. Additionally, the 
value of bandwidth D = 20F is chosen large enough that 
the results in the considered parameter regime do not 
depend on 


III. RESULTS 

Our main focus is on the zeros of the thermopower S, 
which is a rather distinct spectroscopic feature, especially 
when there is a sign change of S. To emphasize the pa¬ 
rameter values where the zeros appear we plot the square 
root of the absolute value of the thermopower and 

use a grayscale scheme where the zero value corresponds 
to white. 


A. Non-interacting case T = 0 


We start by discussin g t he thermopower for the non¬ 
interacting U = 0 case.l^^^*^ If the temperature is the 
lowest energy scale in the system, i. e., kyT F, the 
zeros of the thermopower can be obtained directly from 
the transmission function by using the Mott approxima- 
tionpIHZil which yields the following expression 


Sm 7r\ ^ainr(T) 


(19) 


Then we see that in this limit the zeros of thermopower 
are given by the zeros of the derivative of the trans¬ 
mission, i. e., Sm = 0, when dET{E)\E=o = 0. Using 
Eq. this condition defines the following curves in 

the (F,AT,V^) space: 




sjvi + (ar)2, 


AE 

AE 


1-a^ 




l + a2 

^(1 - a ^) = -(1 + a ^) V , ± + [(1 


(20a) 

(20b) 

a2)r]2. 

(20c) 


From the above expressions we see that for given detun¬ 
ing AT between the levels there can be up to five zeros, 
when the gate voltage Vg is scanned. Three zeros appear 
once the distinct peaks, related to the energies of the sin¬ 
gle particle levels, become resolvable and this happens 
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Fig. 2. (Color online) The thermopower S (left column) and 
conductance G (right column) as a function of the detuning 
AE and gate voltage Vg in the non-interacting case U = 0. 
The asymmetry values are a = 1, 1/2, 1/4, 0 for (a) to (d), 
respectively. For the contour plots the temperature is set to 
/cbT = 0.5F and ± denotes the sign of the thermopower for 
these contour plots. The dashed (blue) curves depict zeros 
of thermopower, given by Eqs. (20), for the case of a small 
temperature F ^ kBT ^ 0. Panel (e) gives a cut of the 
contour plots at the detuning AE = — 8.5F. 


when the separation between the levels is larger than the 
energy scale determined by F. These zero’s curves are 
given by the two solutions ([20^) , and by solution in ([2^) 
with the minus sign for > 0 or the plus sign for <0. 
Additional two zeros appear due to destructive interfer¬ 
ence, and for finite detuning AE ^ 0 this happens only 
when the couplings to the levels are asymmetric, a ^ 1. 
For symmetric couplings, a = 1, destructive interference 


Fig. 3. (Color online) Same as Fig.[^just for the case where 
the tunneling amplitudes have the same sign. The dashed 
(blue) curves depict zeros of thermopower, given by Eqs. (21). 


leads to zero thermopower for AE = 0, but without the 
sign change, because it is a solution of double multiplicity. 
Figure [^ depicts the zeros given by analytical expressions 
Eqs. ( [^ by blue dashed curves for different asymme¬ 
tries a. These lines can be compared to the thermopower 
(in grayscale) where the temperature is compa¬ 
rable to the coupling F. We see that for larger values of 
AE or Vg the zeros, corresponding to single particle lev¬ 
els, for the finite temperature case match Eqs. (20), but 
there are some discrepancies when AE and Vg do not 
exceed Foi* the interference features we see that 

they are sensitive to finite temperature, i. e., depending 
on the coupling asymmetry they appear for larger values 
of AE^ Vg ^ kBT. Fig.|^ shows the calculated data for 
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a fixed AE = —S.5T. For a = 1/2 (dashed blue curve) 
we clearly see the five zeros in thermopower for the finite 
temperature calculation. 

For the case of the same sign tunneling amplitudes 
from Eqs. (18) and (19) we find that there can be up to 
three zeros in the thermopower for T ^ 0 for given AE: 


AE 


AE 

~Y 


1 - a2 


(21a) 

(21b) 


For low energies in the transmission Eq. (18) the destruc¬ 


tive interference corresponds to a zero given by Eq. ( |21b[ ), 
i.e., the gate voltage has to be within the window of en¬ 
ergies given by the detuning, i.e., \Vg\ < \AE\. This 
is in contrast to the canyon of conductance suppression 
situation, where the destructive interference appears for 
\Vg\ > \AE\. Additionally, the canyon of conductance 
suppression exhibits a thermopower zero for \Vg\ < \AE\ 
just because there is a minimum of conductance, and 
there is no necessity for destructive interference for this 
zero to appear. On the other hand, for the same sign 
amplitudes the large values of the thermopower close to 
a zero give a hint for the presence of the destructive in¬ 
terference (see Eig. [^.1^ 

The role of finite temperature can be understood as fol¬ 
lows. In the different limit, when the tunneling coupling 
is the smallest energy scale, i. e., F ^ \E — Ei\^ |E — Ei/ 1 , 
we obtain the following expression for transmission to 
lowest order in F 


T{E) « ttF [S{E - El) + a‘^5{E - Er)] , 


( 22 ) 


which is valid for both Eq. ( |15[ ) and Eq. ( |18| ). Applying 
Eq. ( 8 b) yields the condition for zeros of thermopower 


X cosm y + o?y cosm x = 0 , 
x = Ei/{2kBT), y = Ei,/{2kBT). 


(23) 


The above Eq. ( [M| ) can have only up to three gate voltage 
Vg solutions for given AE, i. e., the features arising due 
to destructive interference for the canyon of conductance 
suppression get washed out due to high temperature [as 
also seen from Eq. ( [ 2 ^] . Additionally, for symmetric 
coupling, a = 1, Eq. (|23|) yields the following condition 
for resolution of two levels. 


e = 


r + 1 
r — 1 ’ 


r = 


AE 

2A)rT’ 


(24) 


which has a solution | AEd ~ S.OQ/ceT.I^Erom Eig.[^ we 
see that both the temperature T and the coupling F con¬ 
tribute to the effective broadening and the two levels be¬ 
come resolvable in the thermopower for AE ^ 2 r -I-S/crT 
in the case of the canyon of conductance suppression. For 
the same sign tunneling amplitudes we obtain that the 
resolution of the thermopower zeros around AE = 0 is 
determined just by the temperature broadening and not 
by F (see Fig. 



Fig. 4. (Color online) The thermopower S and conduc¬ 
tance G evolution as a function of increasing charging energy 
for symmetric coupling configuration a — 1. The charging 
energy values are U = 0, SF, 1 OF for (a) to (c), respec¬ 
tively. For the 2vN thermopower calculations we have set 
/cbAT = O.IF, /cbT = F and for conductance calculations 
we have set feAT = 0, feT = F, F = O.IF. Also a finite 
bandwidth D = 20F is used. The dashed (green) lines de¬ 
note the bare resonances given by Eq. (25). The dashed (red) 
lines show a region of detuning, where the thermopower is not 
defined because of the zero conductance. 


B. Finite 1/ / 0 case 

In order to study the impact of the finite interaction 
U 7 ^ 0 , we first examine a case of equal coupling strengths 
for the two levels, i. e., a = 1. The corresponding ther¬ 
mopower and conductance are shown in Fig.[^ Changing 
the interaction strength U causes a minor change in the 
conductance, while the changes are more visible in the 
thermopower. We note that for 1/ 7 ^ 0 the 2 vN method 
can give negati ve dia gonal reduced matrix elements 
defined in Eq. (A18) and also the linear conductance G 
can become negative, which is unphysical.^^ This is the 
reason that we do not address the interference features 
in the thermopower appearing at large gate voltages Vg. 

The conductance plots in the right column of Eig. 
show, that the conductance peaks for large |AE| appear 
at particular resonances, where the quantum dot ground 
states differing by a single charge cross in energy. These 
are the common sequential tunneling peaks occurring at 


Vg = ± 


U+\AE\ 


(25) 


which are shown by dashed green lines in Eig. Eor 


zero 
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Fig. 5. (Color online) The thermopower S and conductance 
G evolution as a function of increasing charging energy for 
asymmetric coupling configuration a — 1/2. The charging 
energy values are = 0, 2.ST, ST, 7.ST, lOT for (a) to (e), 
respectively. Other parameter values are as in Fig. 


detuning /S.E these sequential tunneling lines are broken 
by a canyon of conductance suppression. 

Now we consider the thermopower in the left column of 
Fig.g Around the sequential tunneling lines there is an 
equal amount of electron and hole tunneling for the given 
resonance and the average energy of tunneling particles 
becomes zeroThis provides a zero in the thermopower. 
As can be seen for large |AF^| there is actually a zero 
close to the dashed green lines, as expected. A further 
zero is occurring in between (at = 0 for the case a = 1 
considered here), where electron and hole tunneling for 
different resonances compensate each other. 

For the charging energy, U < 74 = T + 3/cbT, smaller 
than the effective broadening of the levels, there is a re- 


Fig. 6. (Color online) The thermopower S and conductance 
G evolution as a function of increasing charging energy for 
the same sign tunneling amplitudes Eq. (16) and asymmetric 
coupling configuration a — 1/2. The charging energy values 
are 7/ = 0, 2.5F, 5F, 7.5F, lOF for ^ to (c), respectively. 
Other parameter values are as in Fig.H The (red) horseshoe 
patches denote the region of parameters, where our numerical 
procedure for solving the 2vN equations does not converge. 


gion of detuning, where the level crossing of empty with 
singly occupied and singly with doubly occupied states 
in the dot can not be resolved (see Fig. IT ,b). For the 
corresponding thermopower, there is only one sign shift 
for all gate voltages around zero detuning. For larger de¬ 
tunings AE >Uc — U the levels can be resolved and again 
three sign shifts appear in the thermopower: at the two- 
level crossings and at the electron-hole symmetry point 
Vg = 0. When the charging energy becomes larger than 
the effective broadening there are always three sign shifts 
in the thermopower as can bee seen from Fig. H- 
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-2 0 2 -2 0 2 


/^E/T AE/r 

Fig. 7. The thermopower S (first column) calculated for 
canyon of conductance G suppression (second column) dis¬ 
cussed in Ref. iss The charging energy values are U = 2r 
and = 5r for (a) and (b), respectively. Note that here 
the parametrization tLi = — A/O^tp, tm = \/0Tto tLi' = ^o, 
tjiif = \/0Ato, with to = / (27tef), for tunneling ampli¬ 

tudes is used. Other parameter values are as in Fig. 


The behavior of thermopower and conductance for dif¬ 
ferent U when the coupling is asymmetric a = 1/2 is 
shown in Fig. Here we focus on the thermopower 
zeros, which appear for low values of gate voltage and 
detuning, i.e., we do not address zeros, which appear 
due to interference (as discussed for 7/ = 0 above). The 
canyon of conductance suppression is now tilted due to 
the asymmetric coupling as discussed already in Ref. [54l 
More interestingly, the appearance of the thermopower 
zeros is changed quite drastically: for U < 74 = 5F, 
the curve of zero thermopower becomes strongly tilted 
around /S.E = 0. It takes a characteristic S-shaped form, 
around its center denoted by a box in Figs. [^-c. In¬ 
creasing U pushes two further thermopower zeros (their 
onset is denoted by two asterisks in Figs. [^,c) towards 
the central S-shaped curve. Around critical Uc there is a 
qualitative change (compare Fig. and Fig. B) , where 
the S-shaped curve merges with the other ones. Finally, 
for large U (Fig. |^) the situation resembles the case for 
symmetric coupling, a = 1, addressed in Fig. B- 

The situation when all the tunneling amplitudes have 
the same sign is depicted in Fig. For U > Uc we 
see that the evolution of zeros is qualitatively the same 
as in Fig. However, for the canyon of conductance 
suppression the zero lines have a jump around AE ^ 0 
(see Figs. B and|^). Additionally, for U < Uc the S- 
shaped zero thermopower curve does not appear for the 
same sign configuration. 

Lastly, in Fig. we present the thermopower behavior 
for a system, whose conductance was studied experimen¬ 
tally in Ref. |53l These results may be relevant for future 
thermoelectric experiments on similar type devices, and 
shows the situation when the couplings have more general 
structure than Eq. (13). Note that for given parameters 


in Fig.j^we have the level Ey more strongly coupled than 
the level Ei, which is opposite from the cases studied in 
Figs. [2][5l This results to a mirrored behavior of ther¬ 
mopower and conductance along zero detuning AE = 0. 
Otherwise, the introduced additional asymmetries in the 
couplings between the left and the right lead does not al¬ 
ter the qualitative behavior of the asymmetric coupling 
a = 1/2 thermopower seen in Fig. 


IV. CONCLUSIONS 


The canyon of conductance suppression known from 
Refs. 123 and [54] has been further investigated with ther¬ 
mopower S acting as a probing tool. Zeros of the ther¬ 
mopower are a telling spectroscopic feature, which for the 
system studied in this paper yields information about the 
level coupling asymmetry a (compare Fig. vs. Fig. 
and charging energy U compared to the effective broad¬ 
ening of the levels given by ^ 2F + S/cpF. Additionally, 
for the non-interacting case it was shown that up to five 
zeros can be observed when scanning the gate voltage for 
given detuning. This shows that the thermopower mea¬ 
surement could be useful to resolve features appearing 
due to destructive interference [see Fig. and Eq. (20)]. 
Additionally, by comparing the same sign tunneling am¬ 
plitude thermopower [Eq. (13)] to canyon of conductance 
suppression thermopower it was shown that the canyon of 
conductance suppression has a unique signature around 
detuning AE = 0 (compare Eig. [^vs. Eig. [^. 


zero 
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Appendix A: Second-order von Neumann (2vN) 
approach 

In the 2vN approadP^^^^^^we approximately solve the 
equation 


ih^y=[H,p], (Al) 

by considering the density matrix p elements, which con¬ 
nect the states differing by up to two electron or hole ex¬ 
citations. Such a treatment corresponds to the so-called 
resonant tunneliM approximation in real-time diagram¬ 
matic approach,12Sl and yields an exact current for non¬ 
interacting systems (7/ = 0 for our considered system).!^ 
We note that the 2vN method equations were originally 
derived in Ref. (66] and here we present this derivation for 
convenience with a slightly different notation. We write 
the governing equations in the many-particle eigenbasis 
I a), 1 6),... of the dot Hamiltonian Hy, Expressed 





















in this many-particle basis the tunneling Hamiltonian 
@ becomes 


Ht = 


E 

ab^kl 


{Tba,lW){^ 


a\c^^+\i.c 


)■ 


(A2) 

(A3) 


Tba,i= XI ^*<'(^141®)- 

Here we used the letter eonvention: if more than one 
state enters an equation, then the position of the letter 
in the alphabet follows the particle number (for example 
Nb = Na ^ 1, Nc = Na ^ 2, Na' = Na). In such a 
way the sum restricts to those combinations, where 
Nc = Ni) 1. For our considered system we have four 
many-particle eigenstates 

Eq = 0 , 

El = -V, 


| 0 ), 

|i> = 4lo), 

ii') = 4'io), 

|2)=44|0), E2 = -2Vg, 
and such many-particle tunneling amplitudes 


AE 

AE 


U 
"2 ’ 
U 


(A4) 


Ev =-Vg- — -- 


Tf = 


( 0 Toi/ Toi// 0 

Tio^i 0 0 Ti2/ 

Ti'o/ 0 0 Ti/2,^ 

V 0 ^21/ ^21',^ 0 


/* 


/o ti, 

tu 0 0 

tvi 0 0 -t\^ 

\ 0 tvi -tu 0 


0 

/* 


(A5) 


(A6) 


The density matrix elements are defined as 

P^atbg' = {agWg'), 

where \bg) = \b) 0\g), wit^ 6) denoting the eigenstate of 
the dot Hamiltonian Hd ([2) and \g) denoting the eigen¬ 
state of the lead Hamiltonian ii^LR Here the label 
n provides the number of electron or hole excitations 
needed to transform \g) into \g'). For example, we con¬ 
sider the matrix elements of the type 

Pblb'g = {bg\p\b'g), 

{bg - K\p\ag), 

{dg - K-K'\p\bg), 
pfg-^+^>,b'g = {bg-K + K'\p\b'g). 

Here we have introduced the following notation 

K = k, £: 


rbg — K,ag 
rdg—K—K',bg 


(A7) 


\bg + K) = \b)(S)cl\g), 
\bg-K) = \b)(S)cJg), 
\dg-K-K') = \d)®c^,cjg), 
\bg-K + K') = \h)®c\,Ug). 


(AS) 


(A9) 


By neglecting all the density matrix elements with more 
than two electron or hole excitation n > 2 from Eq. 
we obtain the equations 


+ Ebai,KiPalg+Ki,b'gi~^)^"' 

~ kbg,cig-Ki(~^) 


(AlO) 


Pbg,aig-\-Ki\ -^ai6,«i5 


Na 


^fi^Pci-K,bg = (Ec -E,- S6)pE«,bg 


EEcdi,KiP^dlg-n-ni,bg( 1) 

— f —h 

Pcg—K,cig—Ki\ -^cib,Ki 

_ (_ l)Nairp 

Pcg—K,aig-\-Ki\ -^aiO,/ti5 


Nc 


(All) 


d 

,b' g ^ ^Eb — E^^ E^f — g 

+ EbaunP^llg-K+K'+K,b'g^~E}^"^ 

— f —A/ 

Pbg—K,-\-K',cig — K,\ Ci6 ,k, 

_ ^d] (_l)Nairp ^ 

Pbg—K-\-K',aig-\-K'\ -^aib,K, 5 


(A12) 


^ d2] 


g^Pd^K+K',bg ~ (Ed E^ E^i -E’6)4j-re-re',bs 

+ Tdci,K,p[EK-K'+K,bgi~r^°^ 

+ Tdci,K' P^t!g-K-K'+K' ,bgi~r^‘‘^ 

~ Pdg-K-K.',cig-K(~r ’’Tcib.fc 

(A13) 

Note that all indices with subscript 1 like ai, ci, ki are 
summed over. Additionally, phase factors like (—1)^^ 
appear due to order exchange of the lead operators with 
the dot operators, i.e., c^{\b) (g) |^)) = {—l)^^\b) <^c^\g). 
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Summing Eqs. (AlO) and (All) over all the lead states 
\g) we get 


and when going from Eq. (All) to Eq. (A15) we have 
used 




dt 






T h' 

bci,Ki-^cib' ,Ki 




[ 1 ] 


■ 4 ®"' 


bai,fci 

[ 1 ] 


(A14) 


cb,K ~ (Ec -Efc -E’6 )^c6,k 

+ Tcbi,Kl^[}h-K+Ki + Ecdl,Kl^lilb,-K-Kl 

+ ^\xl,-K+K^Ec-i_b,K^ + , 


(A15) 


[2] [uj I n A "i 

Pbig — K+K\^bg Pb\g — K,+K^bg IV ^ Pb\g—K+K\^bg'> 


- .[0] 


[ 2 ] 


.[2] 


o^-‘ = 0^°' +a-(5 

rcg — K,cig — Ki rcg — K,cig — K ' V-^ ^ / rcg — K,cig — Ki’ 

(A19) 

Here we have also assumed that the electrons in the leads 
are thermally distributed according to the Eermi-Dirac 
distribution / and that this distribution is not affected 
by the coupling to the quantum dots. This assumption 
leads to the following relations for Eq. (All) 


.[2] 


Z^Pbig-K,+K,,bg ~ 

9 

/ ^ rcg — K,cig — K J—i^ cc\i 


(A20) 


{Eb -E,+ E,, - Eb,)^tl,-.+.' 


(Alb) 




{Ea-E^-E^,-Eb)^^2,-.-n' 


(A17) 

where we introduced the following notation 


^[0] _ ^[0] 

^bb' — / .Pba.b'a^ 

9 

^cb,K / j Pca—K,ba\ ’ ^6c, 

9 

$[2] y- [2] 

Ca, — K — K' / j rca — K — K,'Ma'> 


(A18) 


[^SJ 


$[2] 


bb' ^—K+K‘ 


' ^ ^Pbg—K-\-t^'^b'g") 


/« = !ki = (exp[(£'fe - g.i)lh^T(\ + 1) \ 
f—K = 4 fki^ 


for Eq. (A12) 




[ 1 ] 

aib' ,k' ’ 


f-i) 

/ j rcT a — K,-\-K,'— K,',b'a\ / 


Nb 


f 




W/O™ , 

/ ^ rbQ—K-\-K\a-\ q-\-k' \ / 


-f 


and for Eq. ( |A13 ) 

f-E 

/ ^ rcT a — K—K,'-\-K,,ba\ ) 


Nc- 


— f 

^cib,K'J>^^ 


rcig — K — K,'-\-k', bg\ / 


$ 


[ 1 ] 

Cib,K 


f-n 

/taels '—«' 


Nb 


$ 


[1] 

dci ,K, 


d-K, 


(- 1 ) 

b'dg—K—K,',cig—K,'\ ) 


Nb 




(A21) 


/«', (A22) 


For the stationary state we assume the conditions 


d 




9 


db 


d 


[ 2 ] 


d 


db 


which allow to write in terms of as 
db 




Tbai,K^lill)/I k + Tbci,K'^\^ji)/,Kf-K + ^\)ai,K,'^aib' ,k' f k 


bb',-K+K' 


db,—K,—K,' 


Ei^ — E,^/ — Ejy + Ejjf + iTj 

~2dCi,fC^ijb'/fC + ,,/k' — ,^/Tcjb^K/_K + ^lil_^^^Tcj^b,K' f-K 

En + Ei^i — Ed + Eb + ig 


(A23) 


(A24) 
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Here we have added a posit ive in finitesimal = +0 to ensure a proper decay of initial conditions. After inserting the 
above expressions into Eq. ( A15| ) we obtain the integral equations of the 2vN method for the stationary state 


0 = -{E, -E, + Ek + 


+ 


+ 


+ 


+ 


Tcbi,Ki aib,Ki 


^ cdi ,Ki 


En — + E5 + 177 

Tj f -u f T u —Tj f f T u 


Ek. + E^^ — Ed^ + E5 + 177 
f f Tu —Tu f 

cdi,Ki J —Ki^ (llCi,K, ' ^ cbi,K,>^ biCi,K,i cbi ,K J biCi ,K,i ^cdi,K,i 


f-K,Tdi 


Ci,K 


^ ci6,«i 


En — E^^ — Ec E Ec^ + 177 
Tcbi.KxfE ^\}^^^f-KxT]j^ai,Kx—Tcbi,KfK^\)lai,Ki 




^ a\b^K\ 


E,. 


En^ — Ec 


■Ea 


177 


0={E,- E,,Ebl + + Tbeu^EE. 


-<h 


[ 1 ] 


bci,Ki-^cib' ,K,i ^bai,Ki-^aib',K,i- 


[ 1 ] 


(A25a) 


(A25b) 


r 


Additionally, we impose the normalisation condition for 
the diagonal reduced-density matrix elements: 


E*S = i 


(A26) 


The integral equation (A25a) under interest has the 
structure 




(A27) 


It is solved iteratively on an equidistant energy grid E^ 
by having = 2^^ discretization points for our consid¬ 
ered calculations. The zeroth iteration of = E^ is 
determined by making a local approximation, i.e., terms 
of the form which have integrated momentum label 

Ki are neglected. Then the first correction is determined 
as The higher order corrections 

\[i] ’ 


are given by S^'kjn = J2 ki ^'^,'^i^Eln-i *1^® solu¬ 
tion is expressed as -|- make 

iterations up to tt. = 6, which yields good convergence 
for most parameter values. In these iterations we need 
to evaluate a lot Hilbert transform of the form: 


= lj 




-D En — Ei^f ± 177 


D Ek ~ Ek,i 


(A28) 




The principal value integrals are efficiently evaluated on 
equidistant grid with N points using fa st Fo urier trans¬ 
form, which has complexity 0{N\ogN)^^^^ 

Finally, we are interested in the current going from the 
lead £ into the quantum dots, which is given by 




keh 


(A29) 


and which shows that are the energy resolved cur¬ 

rent amplitudes. 


Appendix B: Mapping to generalized Anderson 
model 

In linear response to applied bias V or temperature 
difference AT the Hamiltonian Eq. 0 can be mapped 
to generalized Anderson model by mixing the left lead 
electrons with the right lead electrons, and the dot orbital 
1 with the orbital 1'. The generalized Anderson model 
has the form!^ 

E = E EAa<^ka + F + h-c.) +Un^n^ 

k,a k,a 

+ (s — cr^ cos^ j Tier — -h sin^, 

rr ^ / 


(Bl) 

where rin- = a G {t^i} represents a pseudo-spin 


and h is an effective magne tic fi eld. The Hamiltonian 0 
is expressed in the form of ( |B1[ ) by performing a singular 
value decomposition (SVD) on the tunneling Hamilto¬ 
nian 0 as 




A A ■ 

-L/C ^Rk. 


A 


A = 


thl tLl' 
tm tRif 


= eI 


d, 

di,_ 


0 

0 n 


h.c. 


Rd- 


(B2) 


(B3) 


Using the rotations Ri and the pseudo-spin operators 
are expressed in terms of original basis operators as: 


ffi/c 

= Ri 

Cl/c 


d^ 

= Rd 

di 

7 

Clk 


CRk 


di 


dy 


(B4) 
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For the tunneling amplitudes of the form Eq. (13) we 
obtain 


0 


which corresponds to parallel field configuration. For the 
same sign tunneling amplitudes we obtain 

^-=*(l ~a) = o) ^"’+’ 


R. 


'd,+ 


+1 +n 
Vl + \^—a +1 




- arctan 


2a 


l-a2 


Lastly, for tunneling amplitudes examined in Fig. we 
get: 

^ ^A4) (^0^ 0.54) 


Ri 


( 0.91 

0.41\ 

Rd ~ 

0 

CO 

0 

0.95\ 

1-0.41 

0.91/ ’ 

1 0.95 

0.30/ 


, 0 ^ 0.61. 


In all above cases we have 5 = — V^, h = AF/2. 
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